Strong parity mixing in the FFLO superconductivity in systems with coexisting spin 

and charge fluctuations 
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We study the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state of spin fluctuation mediated pairing, 
and focus on the effect of coexisting charge fluctuations. We find that (i) consecutive transitions 
from singlet pairing to FFLO and further to Sz = 1 triplet pairing can generally take place upon 
increasing the magnetic field when strong charge fiuctuations coexist with spin fiuctuations, and (ii) 
the enhancement of the charge fiuctuations lead to a significant increase of the parity mixing in the 
FFLO state, where the triplet/singlet component ratio in the gap function can be close to unity. 
We propose that such consecutive pairing state transition and strong parity mixing in the FFLO 
state may take place in a quasi-one-dimensional organic superconductor (TMTSF)2X. 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, 
in which the Cooper pairs formed as (fc + Qc h —k + Qc i) 
have a finite center of mass momentum, is one of the 
most fascinating superconducting states. [l|, 0] One of 
the interesting aspects of the FFLO state is the parity 
mixing, i.e., even and odd parity pairings can be mixed. 
Phcnomenological studies have shown that the mixing 
of the singlet and triplet pairings stabilizes the FFLO 
state. |3| Recent microscopic studies have shown that 
the Sz = triplet pairing is mixed with singlet pairing 
in the FFLO state of the Hubbard model on the two- 
leg ladder-type lattice, [Hj and also on the square lattice, 
where d-wave superconductivity is mediated by spin fluc- 
tuations. 0, Q In ref. it has been pointed out that 
the parity mixing stabilizes the FFLO state, even in the 
vicinity of the quantum critical point where the quasi- 
particle lifetime becomes short due to the scattering by 
spin fluctuations. 

Recent experimental indications of possible oc- 
currence of the FFLO state in CeCoIus, Q a 
quasi-two-dimensional (Q2D) organic materials 
such as A-(BETS)2FeCl4 (BETS=bisethylenedithio- 
tetraselenafulvalene) ^ and k-(BEDT-TTF)2Cu(NCS)2 
n3EDT-TTF=bisethylenedithio-tetrathiafulvalene), 
[IQ] and also in a quasi-one-dimensional(QlD) 
one (TMTSF)2C104 (TMTSF=tetramethyl- 

tetraselenafulvalene) [111, have stimulated extensive 
studies in this field. For (TMTSF)2C104 in particular, 
the possibility of the spin triplet pairing has previ- 
ously been sugges ted experimentally for (TMTSF)2X 



iX=PFe, [13|, \1M CIO4 [la 116|), but a more recent 



NMR experiment by Shinagawa et al. has revealed 
that when the magnetic field is low, the pairing occurs in 
the spin-singlet channel, while when the magnetic field 
is high, the pairing state is either an FFLO state or a 
spin-triplet state. Yonezawa et al. have found that 
the onset Tc exhibits a peculiar magnetic field direction 



dependence at high fields, which may be related to 
the occurrence of the FFLO state, where the direction 
of the total momentum of the Cooper pairs can play 
an important role. However, the microscopic origin of 
the transition between high field pairing states like the 
FFLO or the triplet state remains unexplored. 

Theoretically, various studies have investigated the 
possibility of triplet pairing |17l.ll lEHim,!!! and the 
FFLO state [23, [i^l- In particular, three of the present 
authors have previously shown that the triplet /-wave 
pairing can compete with the singlet d-wave pairing in 
the QID system because of the disconnectivity of the 
Fermi surface when 2k f spin and 2kF charge fluctuations 
coexist. 25, 26, [i^] ^kp spin-|-2A;i^ charge fluctuations 



supported from the fact that diffuse X-ray scattering ex- 
periments observe the coexistence of 2kF charge density 
wave(CDW) and the 2kF spin density wave(SDW) in the 
vicinity of the superconducting phase in (TMTSF)2PF6. 
[H, [23] Moreover, we have recently found that this kind 
of triplet pairing due to 2kF spin-|-2A:j;- charge fluctua- 
tions is strongly enhanced by the magnetic field, [s^l 
Then a naive question arises along this line: what hap- 
pens if magnetic field is applied to a system where spin 
singlet pairing dominates at zero field but triplet pairing 
is closely competing ? If the FFLO state emerges, what 
is its nature ? 

Given this background, in the present Letter, we study 
the FFLO state of spin fluctuation mediated supercon- 
ductivity in low dimensional systems, and focus on the 
effect of the charge fluctuations. We find that (i) con- 
secutive transitions from singlet pairing to FFLO and 
further to 6*2 = 1 triplet pairing can generally take 
place upon increasing the magnetic field in the vicinity of 
the SDW-I-CDW coexisting phase, and (ii) the enhance- 
ment of the charge fluctuations leads to a significant in- 
crease of the parity mixing in the FFLO state, where the 
triplet/singlet component ratio in the gap function can 



2 



be close to unity. Based on a calculation on a model 
for (TMTSF)2X, we propose that such consecutive pair- 
ing state transitions and the strong parity mixing in the 
FFLO state may actually be taking place in this material. 

The anisotropic extended Hubbard model [Fig[T] (a)] 
that takes into account the Zeeman effect is given by 



^Vijn.rij. (1) 



Here t 



Uj + hzSgn{a)5ij, where the hopping tij is 
considered only for intrachain (t^) and the interchain (ty) 
nearest neighbors, tx = 1.0 is taken as the energy unit. 
U is the on-site repulsion, and Vij are the off-site re- 
pulsions: Vx-, Vx2, Vxs are nearest, next nearest and 3rd 
nearest neighbor interaction within the chains, and Vy is 
the interchain interaction. We ignore the orbital effect, 
assuming that the magnetic field is applied parallel to 
the conductive x-y plane, thus we assume a sufficiently 
large Maki parameter. 

The bare susceptibilities, bubble-type and ladder-type, 
are written as 
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where ^cr(fc) is the band dispersion that takes into account 
the Zeeman effect measured from the chemical potential 
/i, and /(^) is the Fermi distribution function. 

Within RPA that takes into account the rnagnetic field 
parallel to the spin quantization axis z, [30| the lon- 
gitudinal spin and charge susceptibilities are given as 



Asp 



X 
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-'-T'l where 



X^^) and Xc 



X-^(fc) = [l + xrik)V{k)]xrik)/A{k), (4) 
X^'ik) - -Xo''ik)[U + V{k)]xr{k)/A{k), (5) 
A{k) = [l + x^o^{k)V{k)][l + xrik)V{k)] 

-[U + V{k)fx"o"{k)xrik)- (6) 

The transverse spin susceptibility, in which we ignore the 
off-site repulsions for simplicity, is given as 



(7) 



The pairing interactions from the bubble and ladder 
diagrams are given as 



V,Z{k) = U + V{k) + %-xi;{k) 



[U + 2V{k)Y 



Xch{k), 



VZik) = U^tiik) 



(8) 
(9) 



V,7b{k) = Vik)-2[U + V{k)]Vik)x''^{k) 

-V{k)\^^k) -[U + V{k)f x**(fc),(io) 
VZik) = 0. (11) 

The linearized gap equation for Cooper pairs with the 
total momentum 2Qc (Qc represents the center of mass 
momentum) is given by 

A5:>""'(fc) = ^ Y.lV,Z'ik -q) + VZ'ik + q)] 
q 

x ^(^f;»-//-f-'^7-)V -'(.), (12) 
^a{q+)+^c7'{-q-) 

where q± — q±Qc, (f)'^'^ {k) is the gap function and Aq"^ is 
the eigenvalue of this linearized gap equation. The center 
of mass momentum Q^, which gives the maximum value 
of Aq"^ hes along the x-direction because of the nesting 
of the Fermi surface [HI and Aq^ takes its maximum at 
Qc = (0, 0) because the electrons do not scatter between 
the different directional spins in this pairing channel. 

We define the singlet and the Sz = triplet component 
of the gap function in the opposite spin pairing channel 
as 



cl^ssik) = [q^^^{k)~cj)^Hk)]/2, 

0STo(fc) = [<l>^^{k)+C^^Hk)]/2. 



(13) 



In our calculation, the spin singlet and the spin triplet 
component of the gap function in the FFLO state is es- 
sentially d-wave and /-wave as schematically shown in 
Fig [T] (b) , so we write the singlet {Sz = triplet) com- 
ponent of the FFLO gap <?!>ss(0st°) in Eq. (fTS]) as ^ssd 
{(psTfo), where SSd(ST/") stands for spin singlet d-wave 
(spin triplet /-wave with Sz = 0) pairing. The eigen- 
value of each pairing state is determined as follows. Ag' 
with Qc = (0,0) gives the eigenvalue of the singlet d- 
wave pairing Assd {Sz = triplet /-wave pairing Ast /" ) 
= (0ssd = 0), while A^* with Qc ^ (0,0) gives 
Afflo- Aq^ with Qc = (0,0) gives the eigenvalue for 
the spin triplet /-wave pairing with Sz = +1 {Sz = — 1) 
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FIG. 1: (a) The model adopted in this study. (b) 
The schematic figure of the gap for d-wave(upper) and /- 
wave(lower), where the nodes of gap (blue dashed lines), the 
disconnected Fermi surface in QID lattice (red solid curves). 

First, to make the argument general, we concentrate 
on a simple model with only the on-site C/ = 1.5 and the 
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nearest neighbor repulsion 14 in the a;-direction. When 
Vx is large, 2kp charge fluctuations tends to develop for 
band fillings close to half filling, so we take the band 
filling n — 1.1, where n =number of electrons/number of 
sites. Here we fix the value of ty at 0.5, but this value does 
not have a specific meaning, and qualitatively (although 
not quantitatively) similar results can be obtained for 
other values of ty. The temperature is fixed at T = 0.01 
here. System size is taken as 2048x64 sites here. 

In Figl21 we show the magnetic field dependence of Qcx 
of the FFLO state, the parity mixing rate ^gx /o / (j^ssd and 
the eigenvalues of the gap equation for (a) Vx = and 
(b) Vx = 0.65. Note that we denote the ratio between 
the maximum value of the Sz = triplet component and 
that of the singlet component of the gap function in the 
FFLO state as "0sT/o/0SSd" hereafter. The dominat- 
ing pairing state changes from singlet d-wave to FFLO 
upon increasing the magnetic field for both Vx — and 
Vx — 0.65, but for sufficiently large field, FFLO further 
gives way to the triplet /-wave state with 5'z = 1 for 
Vx = 0.65, i.e., when the charge fluctuations are present. 
The reason why 5*2 = 1 triplet /-wave dominates at high 
fields can be explained as follows. The presence of the 
charge fluctuations suppresses the spin singlet pairing in 



23,l2ll,[2i,[2a[27| 



teraction and enhance the triplet one 
Secondly, an Sz = I triplet pairing state induced by the 
coexistence of spin and charge fluctuations is strongly 
enhanced by the magnetic field applied parallel to the 
spin quantization axis z. [sO] In fact, such a possibility 
of transition from singlet pairing to FFLO, and further 
to triplet pairing upon increasing the magnetic field has 
been phenomcnologically proposed by Shimahara. [18j 
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FIG. 2: (Color online) The /iz-dependence of Qca; (upper 
panels), the Sz = G triplet/singlet ratio in the FFLO 
(/)grpyo/(/>ssd (middle) and the gap equation eigenvalues (lower) 
for (a) Vx = Q and (b) Vx = 0.65. Other parameters are 
U = 1.5, n = 1.1, ty = 0.5, and T = 0.01. 

Let us now look into the nature of the FFLO state. We 
show the Vx dependence of the ratio 0st f° / <^SSd in FigOl 
We see that the mixing ratio increases as T4,, namely. 



the charge fluctuation increases. The strong mixing of 
the triplet pairing component may be expected from the 
fact that the presence of charge fluctuations makes tri plet 
pai ring more competitive against singlet pairing, 0l2l|, 
m, l26l.[27l. [soj This can be seen clearly in inset flgure on 
FiglSl where we plot the eigenvalue of the singlet d-wave 
and Sz — Q triplet /-wave pairing using the formalism 
adopted in ref.[30|, where mixing between the odd and 
even parity pairings is prohibited. 
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FIG. 3: (Color online) The Vx-dependence of the ratio 
0ST !° /4'ssd and the ratio of the eigenvalues between the 
Sz = triplet and the singlet on inset figure. Other pa- 
rameters are the same with Fig[2] 

We now move on to a realistic model for (TMTSF)2X. 
We introduce not only U and Vx, but also other distant 
off-site repulsions Vx2, Vxs and Vy, since the 2kF charge 
fluctuations, which becomes competitive against 2kp spin 
fluctuations when the condition of Vx2 + — U/"^ is 



satisfied, are important in this material. 25|, [2i,|23,[3C 
Here, we fix the repulsions as U ~ 1.7, Vx = 0.9, Vx2 = 
0.45 and Vxs = 0.1, and vary Vy. Other parameters are 
taken as ty = 0.2, T = 0.012, n = 1.5 (3/4 filling), and 
the system size is taken as 1024x64. 

In Fig. m the center of mass momentum Qcx, the par- 
ity mixing ratio (figT f° /4'ssd and the eigenvalues are plot- 
ted as functions of hz in the (a) absence {Vx, Vx2, Vxs and 
= 0) or (b) presence of the off-site repulsions, where 
for the latter case we set Vy = 0.35, for which 2kF charge 
fluctuations are smaller than 2kp spin fluctuations. As in 
the previous case for the simple model, FFLO dominates 
over singlet d-wave pairing upon increasing the magnetic 
fleld, and the FFLO state further gives way to the 5*2 = 1 
triplet /-wave state in the presence of the off-site repul- 
sions. As shown in the middle panels in Fig [41 a strong 
mixing of singlet and Sz = triplet pairing components 
takes place in the FFLO state, especially when the charge 
fluctuations are strong. 

Finally, we show in Fig. [5] a phase diagram for the 
pairing competition in the off-site repulsion Vy versus 
magnetic field hz space obtained by comparing the eigen- 
values of the gap equation for each pairing channel. 
The size of the symbols denotes the magnitude of the 
eigenvalues, and "SDW-fCDW" means that both spin 
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FIG. 4: (Color online) The /i^-dependence of Qca: (upper), the 
triplet/singlet ratio (/)gx/o/(/>ssd (center) and the each eigen- 
values (lower) in (a) the absence and (b) the presence of the 
off-site repulsions. 



and charge susceptibilities have divergently large values. 
The phase diagram shows that for large enough Vy, the 
singlet^FFLO(with strong parity mixing) ^triplet tran- 
sition takes place upon increasing the magnetic field. 
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FIG. 5: (Color online) A pairing phase diagram in Vy-hz 
space, where the green circles represent the singlet d-wave, 
the red squares the FFLO, the blue triangles the = 1 
triplet /-wave and the black filled diamonds SDW-I-CDW. In 
the shaded area, the eigenvalue is represented by the contours 
and the hatched line area is the SDW-fCDW instability. 

To conclude, we find that (i) consecutive transitions 
from singlet pairing to FFLO and further to 5*2 = 1 
triplet pairing can generally take place upon increasing 
the magnetic field in the vicinity of the SDW+CDW co- 
existing phase, and (ii) the enhancement of the charge 
fluctuations leads to a significant increase of the parity 
mixing in the FFLO state, where the triplet/singlet com- 
ponent ratio in the gap function can be close to unity. 

We raise (TMTSF)2X as a candidate material for such 



consecutive pairing state transition and strong parity 
mixing in the FFLO state to take place. In fact, as men- 
tioned in the introductory part, the experiments suggest 
the presence of low field and high field pairing states. 



where the former occurs in the spin-singlet channel. 11 1 
As for the high magnetic field pairing state, Yonezawa et 
al have shown that for a magnetic field parallel to the a 
axis, only the clean sample exhibits an upturn of the T,. 
curve in the high magnetic field regime above 4T, which 
suggests the presence of a pairing state sensitive to the 
impurity content. I4I Between 4T and the Pauli limit 
of around 2.5T, there seems to be a different high field 
pairing state, in which superconductivity is stable against 
the impurities, but is very sensitive to the tilt of the mag- 
netic field out of the a-h plane. The bottom line of these 
experiments is that there may be three kinds of pairing 
states, i.e., one low field state, and two high field states. 
The correspondence between these experimental obser- 
vations and the present study is not clear at the present 
stage, but the appearance of three kinds of pairing states 
is indeed intriguing. It would be interesting to further 
investigate experimentally the possibility and the nature 
of two kinds of high field pairing states. 
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